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We investigate quantum fluctuations on a de Sitter (dS) brane, which has its own thickness, 
in order to examine whether or not the finite thickness of the brane can act as a natural cut- 
off for the Kaluza-Klein (KK) spectrum. We calculate the amplitude of the KK modes and the 
bound state by using the zeta function method after a dimensional reduction. We show that the 
KK amplitude is finite for a given brane thickness and in the thin wall limit the standard surface 
. divergent behavior is recovered. The strength of the divergence in the thin wall limit depends on 

' the number of dimensions, e.g., logarithmic on a two dimensional brane and quadratic on a four 

' dimensional brane. We also find that the amplitude of the bound state mode and KK modes depends 

on the choice of renormalization scale; and for fixed renormalization scales the bound state mode is 
(— I . insensitive to the brane thickness both for two and four-dimensional dS branes. 
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I. INTRODUCTION 
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00 ' Recent progress in string theory suggests that our universe is not four-dimensional, in reality, but is a four- 
' dimensional submanifold, called a "6rane" embedded into a higher dimensional spacetime, called the ^^bulk" . The 
1^ model which was proposed by Randall and Sundrum (RS) succeeds in the localization of gravity around the brane 
' due to a fascinating process, i.e., through the warping of the extra-dimension Q]. This model has been given phe- 
<^ , nomenological grounds from various aspects of higher-dimensional theories of gravity, e.g., in terms of cosmology, 
^ ' higher-dimensional black holes and the AdS/CFT correspondence (see e.g., j2ni3 for reviews). 

In spite of much effort by various authors, there seem to be few theoretical predictions of such a braneworld with 
^ •■ cosmological observations. In the RS type brane model, information from a higher-dimensional gravity theory is 
carried by Kaluza-Klein (KK) modes, which are massive modes from the observers' point of view, who are living on 
^ the brane. Thus, if KK modes are detected in future observations this may uniquely determine whether or not our 
_ universe is a brane. 

, The quantum fluctuations that are usually assumed to be produced during the inflationary phase are considered as 
d ' natural seeds for the perturbations, which form the various cosmological structures and produce the CMB anisotropy 
we see today. Thus, in the braneworld scenario, the quantum fluctuations may imprint information from the KK 
modes on today's cosmological observations. 

In some braneworld models, the inflaton whose dynamics induces inflation on the brane is naturally set into the bulk 
as a result of a dimensional reduction of some higher-dimensional gravitational theory, e.g., dilaton or moduli fields, 
as discussed in In various bulk inflaton models, quantum fluctuations have been well examined jl, [6|, H, ■ 

However, it is well-known that a Casimir surface divergence remains on the brane even after UV regularization. 
This type of divergence prevents us from evaluating quantum fluctuations exactly on the brane, though there are 



approaches to deal with such a problem, e.g., see |8|,P|. In all previous cases mentioned above it was assumed that 
the brane was inflnitesimally thin, which is inherited from the original proposal by RS. 

From a more realistic point of view; however, it is rather natural for the brane to have a thickness. Furthermore, 
string theory, which brane world is based upon, has a minimum length scale, i.e., the string length scale Ig — \/a', 
where a' is the inverse of string tension. This supports the possibility that, in reality, the brane has a flnite thickness, 
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rather than being infinitesimally thin. Various thick brane models have been discussed, see e.g., [Tol Hit ITsl IT^ 
Ellllllllillla Then, we are interested in whether or not "t/ie finite thickness of the brane, if it exists, acts as a 
natural cut-off for the KK spectrum." In this article, we discuss this speculative idea by investigating the behavior of 
the KK modes in an explicit thick de Sitter (dS) brane model. 

We consider a simple thick dS brane model which is supported by a bulk scalar field (p with an axion-like potential, 
originally discussed by Wang [23, E| ■ Then, wc introduce another scalar field x as a bulk inflaton 0, Q type field and 
investigate its quantum fluctuations. We use the zeta function method in conjunction with the dimensional reduction 
approach developed in j^]. We show that the brane thickness regularizes the UV behavior of the KK modes and 
compare its amplitude with that of the bound state mode, which is relevant from the observational point of view. 

This article is organized as follows: in Sec. II, we introduce a thick dS brane model supported by a bulk scalar 
field. In Sec. Ill, bearing the bulk inflaton model in mind, we introduce another quantized scalar field and make the 
necessary preparations to evaluate the amplitude. We take the dimensional reduction approach and calculate the 
amplitude using generalized zeta functions. In Sec. IV, focusing on a two dimensional dS brane (three-dimensional 
bulk), we calculate the amplitude of quantum fluctuations and show that the brane thickness acts as a natural cut-off 
for the KK spectrum. The surface divergence in the thin wall limit is logarithmic. We also show that the bound state 
amplitude depends on the renormalization scale; however, for a fixed renormalization scale it is almost independent of 
the brane thickness. Then in Sec. V, we discuss quantum fluctuations on a four dimensional dS brane (five-dimensional 
bulk) and obtain results similar to that for the two-dimensional case. The leading order surface divergence in the 
thin wall limit is quadratic in this case. In Sec. VI, we summarize our results and mention future work related to 
these issues. In Appendix A, we derive the fluctuation amplitude of the bound state on two and four dimensional 
dS branes, respectively. In Appendix B, we demonstrate the classical stability of the thick brane model both against 
tensor and scalar perturbations. 



II. A THICK DE SITTER BRANE MODEL 



We consider the Einstein theory coupled to a bulk scalar field, 

S = - J d'+'xV^^ R -(50)'-2y(0)) 



where the potential of the scalar field is given by the axion like form |20l 121 

2(1-<t) 



Vi(b) = cos 



^90 



(2.1) 



(2.2) 



Note that we shall set k^,^^ = 1 in this article. 

We shall assume a static configuration, namely (j) depends on only the bulk coordinate and make the following 
metric ansatz 



ds' = b{zy{dz^ + -f^^dxf'dx'') , 



(2.3) 



where 7^^ denotes the metric of d-dimensional de Sitter (dS) spacetime. Following the above ansatz, we obtain the 
Einstein equations 



2 '\b, 
-b" 1 
T+2' 



2 2 



b^V 



and the field equation for the scalar field is 



(2.4) 



(2.5) 



where the prime ' denotes the derivative with respect to z. Note that only two of these three equations are indepen- 
dent. For this potential, we find the solutions 



b{z) — ( cosh ( 



So sm 
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tanh 



(2.6) 
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where 

^^- (.-DtiTld-i).] ' ^o = V(^-iMi-.). (2.7) 

This solution represents a dS domain wall whose energy is localized at z = 0, i.e., the center of the wall. The parameter 
CT has the meaning of the thickness of the wall (brane) from the physical point of view. In order to keep the positivity 
in the square root, we should restrict the range to 

< (T < 1 . (2.8) 



III. QUANTIZED SCALAR FIELD PERTURBATIONS 

Our purpose is to discuss the quantized scalar field perturbations on a thick, inflating brane model. We achieve this 
by introducing another scalar field Xi which is coupled to the domain wall configuration and its fluctuations. Hence, 
we add the action of the scalar field x to the original action Eq. (|2.I|I . i.e., 

S = \j d''+'x^g(T -{d^y - 2V{^)) + \ J d'+'xV^(-{dxY e x') , (3.1) 

where ^ is the scalar curvature coupling. 

As discussed in [23, 01 the coupling of the field x to can be ignored because its backreaction to the domain 
wall geometry is only important at higher order, 0{x^)-^ This assumption allows us to treat the x-field contribution 
perturbatively. The minimally coupled case, ^ = 0, will be of particular interest, because the perturbation equations 
are very similar to those for tensor perturbations of the metric (see Appendix B). 



A. Dimensional reduction approach 



We shall evaluate the amplitude of the quantum field x based on a dimensional reduction of the higher dimensional 
canonically quantized field. This method has been already discussed in Q and we refer the reader to this reference 
for more details. 

In this method, the action of x is rewritten as 

Sx^lJ d''+\^gx{^d+i-^^''R^)x, (3.2) 

where we set a regulator boundary at z — L in order to obtain a well-posed quantum field theory on the dimensionally 
reduced spacetime. Then, the bulk modes become discrete and the solution is written as 

x{z,xn=Y.P-('^)'fir.{xnH'/\ (3.3) 

n 

where ipn has the dimension of a scalar field in the d-dimensional dS spacetime. Due to the maximal symmetry of dS 
spacetime, we can integrate out the dependence on the transverse directions, x^, assuming that the vacuum respects 
the dS invariance. Hence, we shall drop it in the amplitude. 

Integrating the action with respect to z, it is reduced to the summation of theories of a c?-dimensional massive scalar 
field with mass m„: 




n 



where we employed the normalization condition 

2/ {Hdz)b''-\z)FgJz)Fq,Jz)^Snn' . (3.5) 

"'0 



^ These works use the methods developed in l25l . 
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FIG. 1: The potential for a minimally coupled bulk scalar 
field perturbation is shown as a function of Hz for a four- 
dimensional dS wall. The thick, thick-dashed and dashed 
curves correspond to the cases of cr = 0.01,0.05,0.1, re- 
spectively. 



FIG. 2: v/a IS plotted as a function of thickness, a. The 
thick, thick-dashed and dashed curves correspond to the 
cases for ^ = 0, 3/32, 3/20, respectively. 



Note that the multiplying factor of two is due to the Z2-symmetry. The mass-squared m„ is given by 

id - 1)^ 

ml = qlH' + ^^H' . (3.6) 
We introduce a new function fq^{z) := b^'^^^^^^{z)Fq^{z), which obeys the Schrodinger like equation 

- (z) + V{z)H^U^ = qlH^U. {z) , (3.7) 

where 

V{z) -(e. - 0{d{d - 1) + ^} . (3.8) 

cr J cosh [Hz /a) 

For the minimally coupled case, ^ = 0, this potential reduces to the one for the tensor metric perturbations (i.e., gravi- 
tons), Eq. (|B3I) . The reason that we choose a massless scalar field as a probe in this article is to obtain some insight 
into the graviton case (and also for technical simplicity). In Fig. 1, we plot the potential for the d = 4 case explicitly 
for ^ = 0. It is evident that the potential becomes deeper for smaller values of a. 
The solution of the KK modes is 

/,„ (z) = CiPr- {^) + ^2^-'""" [x) , (3.9) 
where Plf{x) denotes the Legendre functions of the first kind, x :— tanh(iJz/o-) and 



1/ := 



+ 4(ec - Oid{d - 1)(t2 + 2d<j) - 1 



(3.10) 



The coupling 



= ^ (3.11) 

denotes the conformal coupling strength, e.g., for the d — 2 case = 1/8 and for the d — 4 case = 3/16. In this 
article, we restrict the coupling to the range < ^ < ^c- 
The mass of the bound state mode is given by 

iv 

go = - , (3.12) 

(7 

which has a maximal value of {d— l)/2 at ^ = 0. For ^ < 0, the bound state becomes tachyonic and non-normalizable. 
For ^ = 0, it is the zero mode and for C = it becomes the lowest mass KK mode, irrespective of the choice of a. In 
Fig. 2, we plot v/a as a function of a for several choice of ^. We find that this ratio is almost independent of a. Note 
that the ^ = case is equivalent to the bulk tensor perturbations. 
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B. The zeta function method 



Given the functions fq^{z), the vacuum expectation value is defined by 

- P^E/lW«(^))' (3-13) 

where the factor of two is due to the Z2-symnietry. From now on, we shall discuss the quantized field theory in 
Euclidean space, i.e., the metric is 

ds^ = b'^{z){dz^ + dnl) , (3.14) 
where dfl^ is the line element of S"^ with unit radius, whose volume is given by 

2 ~ 

Thus, in order to consider the quantum fluctuations of a d-dimensional field, we assume that the vacuum is given 
by the Euclidean vacuum, which corresponds to the dS invariant, Bunch-Davis vacuum in the original Lorentzian 
spacetime. 

For the d-sphere, S'^, any local quantities are related to global ones by simply dividing by the volume of the sphere (a 
property of maximally symmetric spaces; see 0). Thus , we are particularly interested in the local vacuum expectation 
value as only a function of z (one non-trivial dimension), implying 

Kn{t)^ [ d'^x^Kr,{x,x;t), (3.16) 



where Kn is the dS heat kernel for each mode n, see Thus, due to the maximal symmetry of dS space, the global 
heat kernel is simply related to the local one by 

Kn{t)='^K„ix,x;t) . (3.17) 
At this stage it is convenient to rescale the amplitude as 

{x'iz))^ ^%l^'U x'iz)), (3.18) 

where overall factors can be restored at the end of the calculation. Now we may sum up all the KK modes in Eq. 
(|3.13l) : however, as is well known, a naive summation over all the KK modes gives rise to unwanted divergences. 

To deal with such a problem, we construct the local zeta function, C(z, s), along the lines of reference 8], where the 
parameter s is initially assumed to be s > {d+ l)/2 in {d+ l)-dimensions. Once we obtain such a zeta function, after 
analytic continuation to s — s- 1, we end up with 

{f{z))^\imCiz,s), (3.19) 



C>:*) ^475TT^C(^'^) = -FTT / dtt^-'K{z,t). (3.20) 



where 

K{z,t) is the local heat kernel defined as 

oo 

K{z,t) = 2Y,fgMKn{t), (3.21) 

n=l 

and 

oo 

Kn{t) = Y,d, e-[9i^+0-+i/2)^]ff^t ^ (3.22) 
where dj is the degeneracy for each mode j given by 

{3+d-2) \ 
j\{d-\)\ 



d,^i2j+d-l)'-^LZ__±^ (3.23) 



is the global heat kernel for each KK mode. Note that the dimension of ( is slightly different to the case discussed in 
0, because of a difference in dimension of the warp factor. 
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C. Contour integral representation of the local zeta function 

First, as a resolution to the subtle nature of the continuous modes, we introduce another boundary at z = L. This 
then enables us to evaluate the zeta function using the residue theorem, based on certain assumptions relating to the 
zeros of the function in the contour. Then, after constructing such a zeta function we show that we can take the one 
brane limit L ^ oo in a well defined manner.^ 

The solution for the scalar field perturbations in general dimensions is given by 

fqiz) = Nq(^aqPnx) - PqRTix)) , (3.24) 
where for convenience, we choose the second solution R]^" [x) to satisfy 



1 — x^ 



where x — t'dnh{H z / a) . There are several candidates for R]^'^ {x) such as 

£H9^±ii±i)Qj,«^(,), ^ " , Pr^-(x), (3.26) 
T{iqa + v+l) ^ ^ ' 2i sinh(7r(7cr) ^ ' ' ^ ^ 

and so on. For now we do not need to specify the explicit form of the second solution R]^'^{x), but only use the 
property of the Wronskian in Eq. H3.25(l . 

To be specific, let us consider the case of Neumann boundary conditions. The boundary conditions at the center of 
the thick brane and the second boundary are respectively 

/^(z)U=o-0, f'^{z)U^L^O. (3.27) 

Note, the thick brane is not a boundary, we just fix the z derivative of the field at a point to obtain a well-posed 
eigenvalue equation. From these boundary conditions, we get an equation which determines the KK mass spectrum 
as 

Pr'{0)RT'{xL) - Pr'{xL)RT'{0) = . (3.28) 
We denote the solutions for the eigen-equation as (n = 1, 2, 3, • • •) whose eigenfunctions are 

/,„(z) = 7V,„ (a,„P;«"-(a;) - /3,„i?^^(a:)) , (3.29) 

where 

Rl^"'"{0) RJ-^'ixL) 



(3.30) 



We assume gi < 92 < 93 < • • respectively. Note that the final equality is satisfied only for q = qn- Without loss of 
generality, we can choose — R]^'^'{Q) and (3^ — Pl'^'"{Q). We shall also require the normalization constant for n-th 
mode which is found to be 



7V-„2 ^ 2j\H dz) (a,„ Pl'i-^ix) - /3,„ i?^- {x)) ' 



1 i?j?"'"'(0) 



d. 



^jRT'io)pr'ixL)~pr'imT'ixL)) . (3.31) 

^ / Q=Qn 



Note that in the second step we used Eq. (3.8) and integrated by parts. 

Now we have all the necessary tools to calculate the zeta function by applying the residue theorem as follows: from 
the equations given above, the normalized mode functions can be written as 



The same approach cannot be used for the one-loop effective action because it is a global quantity, e.g., see the discussion in |^. 



7 



where, 

F{q) = -(^RT'{0)Pr'{^L) - Pr'{0)RT'{^L)) (3.33) 

and 

G{q,z) = (i?r'(o)^r(^) - pr'mTi^^)) (i?r'(^L)^r(^) - pr'i^Lwnx)) . (3.34) 

This form is essential in order to apply the residue theorem. Whence, the zeta function can be written as a contour 
integral in the complex u plane 

°° °° d P (z) 

_ 0,,2(s-l) "'JJq^y^) 



00 oc 



2 2{s-i) crqnGiqn, z) dj 

2(8-1) I du auG{u,z) ^ dj_ 



- fc 2m F(u)ul,„ [«2 + + ^)2].^2. ' (3.35) 

where the poles at u — qn are on the positive side of real axis and therefore, the closed contour C has to be taken 
around the positive real axis in general. Note that we have introduced a mass scale fi to keep the dimension. This 
term is in fact the renormalization scale and groups with any divergent terms in the expression for the amplitude. 
Then, given the fact that there are no poles in the complex u-plane, besides those on the real axis, we can naturally 
deform the contour C into C (see Fig. 3) 

C(z a)-2u'^^-'^ i ouGju^z) ^ d^ 

C(z, .) - 2m ^^^^1^^^^ ^ ^ ^^)2YH^s ' (3.36) 

which is composed of a line parallel to the imaginary axis with a small real part and a large semi circle on the positive 
real half of the complex plane, which is depicted in Fig 3. As we mentioned previously, initially keeping s larger than 
(d + l)/2, the contribution from the larger semi-circle becomes negligible. 

A similar approach has been used e.g., in 26j for infinitely thin Minkowski branes in a bulk AdS space; however, 
in our case the contour we have to construct is complicated by the presence of the poles which come from the bound 
state; and as we shall see, it will be convenient to evaluate the bound state contribution separately. Therefore, as it 
turns out, we shall only focus on the total amplitude from now on.'^ 

In particular, we are primarily interested in calculating the mode functions on the brane at z = (x = 0), i.e.. 



G{u,z) 



F{u) 



z=0 



(i?™'"(0)-P,5'"'(0) - P™'"'(0)i?™''(0))(i?t"'"(a;L)^T''(0) - Pr'"(2;L)i?r''(0)) 

-(i?jr'"'(0)P™'"'(xL) - PI'^'^'{0)R^'"{xl)) 
P™'"'(xi)i?™'"(0) - Ri'"'{xLyP^'"'iO) 



{R])""{Q)PI'^'"{xl) - Pi,'""{0)RlJ"''{xL)) 

Pr'{xL)P-''"'{0) - P-'""^(a:L)Pr"(0) 
(P,-™"'(0)P»-'(a;L) - Pr"'(0)P."™"'(a;L)) ' 



(3.37) 



where in the first step we used the Wronskian relation Eq. H3.25II and in the final step we specified the second mode 
function as 

2% smh(7rgcr) 

Two types of decomposition are possible: 

Pt'""{xL)P-''"'{Q) - P-''""{xl)PT''{Q) 
P.-'""'(0)P™-'(xl) - P™'^'(0)P,-™"'(xl) 



^ The zeta function for the zero mode is discussed in Appendix A for d = 2 and d = A. 
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FIG. 3: The contour C, used to evaluate the KK am- FIG. 4: The contour C, used to evaluate the total am- 
plitude. The poles on the real axis qi (i = 1, 2, ■ ■ •) cor- plitude. 
respond to the KK modes, while those on the imaginary 
axis correspond to the bound state. We can deform C 
into C' because there are no poles in the complex plane 
besides those on the real and imaginary axes. The closed 
contour depicted by the dotted line is used to evaluate 
the bound state amplitude. 



P^"'^(0) 2isinh(7riicr) 1 

pm<T/(o) ~ P»'^'(0) ^ P™'^'(0)P^-™'"'(xl) - P™'^'(a;L)P^-™'"'(0) 

It is important to note that the second term on the second line is neghgible in the a;^ — *■ 1 limit on the upper half 
of the complex u-plane, while the second term on the third line is negligible in the same limit on the lower-half of 
complex u-plane. Thus, in the single brane limit we use the first term on the second and third lines as the single 
brane propagator on the upper and lower half of the complex plane, respectively. 

In the single brane propagator given above, P^"'^(0)/P^"'^'(0) has poles that are situated on the negative imaginary 
axis, corresponding to purely decaying modes, plus the bound state contribution at u = iv/a. However, as we 
mentioned above, P^"'^(0)/P™'^'(0) is used for the upper half of the complex plane and we need not worry about the 
purely decaying modes. Thus, we only need to deal with the bound state mode at m = iv /a in the calculation of the 
KK amplitude. Similarly, the exact opposite occurs for Pj7™°"(0)/P^™°"'(0) and we only need to deal with the pole 
at It = —ivju. 

The remaining problem then concerns the avoidance of the bound state poles at u = ±iv / a. We avoid the bound 
state poles by deforming the contour to C", as depicted in Fig. 3, when we evaluate the KK amplitude. However, this 
contour gives a non-zero contribution (from the bound state poles) when taking the Cauchy principal value on the 
imaginary axis. This contribution simply corresponds to the subtraction of the bound state from the total amplitude; 
we can calculate the bound state amplitude separately, see Appendix A. Thus, it will be rather convenient for us 
to shift the contour over the upper pole to C, as depicted in Fig. 4. This is equivalent to adding the bound state 
contribution with a counter-clockwise contour (the closed dotted line in Fig. 3) to C". Then, by integrating along 
the contour C and subtracting the bound state contribution, we can obtain the desired KK amplitude. This is the 
approach we shall take to evaluate the KK amplitude in this article. 



(3.39) 
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IV. KALUZA-KLEIN AMPLITUDE: d = 2 CASE 

To demonstrate the method discussed in the previous section as simply as possible, we first evaluate the amplitude 
of the quantum fluctuations on the brane for d — 2. That is, we construct the zeta function for the case of the 
two-sphere in the transverse dimensions with one non-trivial bulk dimension. 



A. Amplitude of the KK modes 

The zeta function for total amplitude at the center of the wall is 

r du cyuG(u,Q) ^ (j 1/2) 



C(0,s) 



4^^ 



27ri F{u) 



E 



u2 + {j + 1/2)2]'* i72s 



4^- 



2{s~l) 



■sin[7r(s- 1)] P 



3=0 



E 



(j + 1/2) 



j^(eW2[/) ^ [c/2 _ (j- + 1/2)2 



(4.1) 



where U = e ^^/'^u and we use the property 

G(e'^^/2C/,0) _ G(e-"/2C/,0) 



(0) 



i^(e"/2[/) F(e-W2f7) 2P^'^'"{0) 



(4.2) 



Here, Roman "P" (not to be confused with the Legendre function of the first kind) means taking the Cauchy principal 
value in order to deal with the pole at U — v j a. In Fig. 4, the contribution from the anti-clockwise semi-circle around 
u = iv/a cancels with that from the clockwise semi-circle around u = —ivju. 

In the following, we shall divide the integral into two; i.e., for U > ^vja which we denote as the "UV piece" and 
that for < [/ < 2i^/(T which we denote as the "IR piece." We emphasize that the reason for this splitting is solely 
for technical reasons and that the choice of division has no physical significance. We can set the split at any value of 

0(l). _ ^ ^ ... 

To begin with, for the UV piece we will use the following asymptotic expansion formula, e.g., see Ref. |27l |. for large 
[/, i.e., 



E 



- [[/2-(j + 1/2)2]^ 2 



1 



1 



1 



r(.) 



E 



r(j 



1) 



c/-'^aaCH(-2j, 



Then, for the IR piece we employ the standard binomial expansions: 



a=l/2 



E' 



3 + 1/2 



E^^fe^t^''Ci/(2s-t-2J-l,^' 



^t^J(, + 1/2)2 -C/2]' 

which is valid for the range < C/ < 1/2; while for the range 1/2 < [/ < ^vjo < 1 we must use pvl l2g 

-2s-2,/N 



Ett- 



3 + 1/2 



1 



J!r(s) 



f;£ii±i)c/2./(2C^(2. + 2J-l,i) 



(4.3) 



(4.4) 



U [(^' + V2)2 - c/2] ^ 2\{\-u^Y • j^^ J!r(s) 

Then, the total amplitude on the center at the wall is given by the summation of both pieces 

C(0,s) = Cuv(0,s) + Cir(0,s). 
First, let us consider the analytic continuation of the UV piece 
2/^2(^-1) sin[7r(s - 1)] 



(4.5) 



(4.6) 



Cuv(0,s) 



7riJ2^ 

^^aG(eW2t/o)^_,^^3 



2l//(T 



F(eW2[/) 



a=l/2 



(4.7) 
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Given the following relation |29j 

p-^-(o) _ 1 r(-| + ^^)r(| + i^ + i) ^^^^^ 



p-^"'(o) 2r(| + i^ + i)r(-f + ^ + 

and by employing the asymptotic expansion for large U of the Gamma functions ,29| we find the following asymptotic 
series, which in d-dimensions is 



-T.^^i'^^OU-'-'', (4.9) 



where is given by Eq. H3.10|l and 



(2 + aid - l ))(l + d(-l + 4g)) ^ 

80-2 2(7 
(2 + aid - 1))(1 + d(-l + 40) (8 + 6a (1 + d (-1 + 4^) + 3a^ (-1 + d) (1 + d (-1 + 4^)) 



= ^ _L^(o), (4.10) 



128ct4 



The subtraction of the wq term just corresponds to that of the trivial background, whereas the wi term corresponds 
to the tadpole graph, see |2^ l24l|. Here, we require only the term wo, in order to regularize the d = 2 case. For the 
d — A case, terms up to wi are required.'' 

Thus, after analytic continuation to s — > 1, we obtain the UV amplitude 

F-^^-(O) woia,i)^ , _j2v 



p-^-'(O) u 



As for the IR piece it is already finite in the limit s ^ 1; however, because of the poles on the imaginary axis we 
make the principle value prescription, i.e., 

cik(o,.) = sin[.(.-i)]P r d^-^^/-";;^' v ^^-•-c.(2. + 2 j - 1, i) 



-1) pAU/<J 

n , — ^ sin[7r(s — 1)1 P / dU , , 



2S-2JN 

where if 2i'/ct < 1/2 the second term is to be dropped. Then, given the Laurent expansion of the Hurwitz zeta 
function 

C//(2s - 1, ^) = ^(^^ - ^(1/2) + Ois - 1) , (4.13) 

we find that there is only a contribution from J = in both terms. Thus, in the limit s — > 1, we obtain 

7G(e^''/2t/,( 
F(e'^V2[/) 



^"'/'^ dU f Ti-v/2 + t/g/2 + l)r(;//2 + Ua/2 + 1/2) _ 1 v Tiv + \l2) \ 
U -v/a\ Tiv/2 + [/cr/2 + l)r(-i^/2 + Ua/2 + 1/2) ^ ^a r(z^ + 1) ) 



(4.14) 



* In practice, for better numerical convergence we subtract off more terms than are required to regularize the theory; thus, we include wi 
for d = 2 and W2 for d = 4, respectively. 
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where in the final step, we used the fact that 

Xo Jo X~ Xq Jq X~ Xq 

'''\xJ^- Tdxl^, (4.15) 

X- Xo Jo X- Xo 

where f{x) is an arbitrary regular function. The second term, which is equal to zero, eliminates the singularity at 
x = Xo in the first term. This technique will also be used for the d = A case. 
Finally, we obtain the regularized total amplitude 

(x'(0)>tot = Irni (Cuv(0,s)+Cir(0,s)) . (4.16) 

As discussed in the preceding section the KK amplitude is obtained by subtracting the bound state amplitude, which 
is evaluated in Appendix A, 

(x'(0))kk = (x'(0))tot - (x'(0))bs . (4.17) 

Interestingly, the total amplitude (x^(0))tot does not depend on the renormalization scale ^, whereas as shown in 
Appendix A the bound state contribution (x^(0))bs does depend on it. Thus, the KK amplitude (x^(0))kk will also 
depend on ^ as can be readily seen from Eq. (|4.17|) . In other words, the dependence on ^ in the bound state and KK 
contribution cancels when they are summed up. 



B. Results of numerical calculations 

The total amplitude (x^(0)) tot is shown in Fig. 5. For small thicknesses, the UV piece dominates the total amplitude. 
The leading order divergent behavior can be estimated as follows: by changing variables from U io x = Ua, the UV 
piece can be written as 

ffMin,<„v(0„s)^-2(i£.„(|;| + i) +5:) , («S, 

As shown previously (see Fig. 2), ly/cr is almost independent of a and ly = 0{a) for tr <C 1. Then by Taylor expanding 
the Gamma functions in Eq. (|4.8|l about ly we find that 

P,7^(0) 1 V f ,x^ ,x + l^ ,x .\ 

^ - -- + ^ 2 ) - m-^) + ^( 2 + 1)) + O(a^) . (4.19) 

Therefore, 

lim Cuv(0, s) = -^1^^ dx (^^(|) - 2^(^) + VXf + 1)) + 0(a"). (4.20) 

In the case of d = 2 the divergence arises only from the leading order. Furthermore, for a; 3> 1 the integrand behaves 
as x^^ and thus, the contribution from the upper bound vanishes. However, in the opposite limit, a; <C 1, 

^(f ) - + + 1) - -- + ("27 - 2V;(l/2)) + 0{x) , (4.21) 

2 2 2 x 

where 7 = 0.57721 • • • is Euler's constant, and therefore 

lim Cuv(0, s) = - — In (2i/) + 0{a°). (4.22) 

s— >1 (7 

Thus, we find a positive logarithmic divergence in the thin wall limit. 

The amplitude of the bound state is derived separately in Appendix A. Here, we recapitulate the final result, 

1 / 1 \ 2 7 1 / 1 ^ 

H^ix'm^s = -(21n(|:)-2V;(l/2)-^e,o(^2j + E ' Ci^(2 J + 1, 2 ) " ^0 2 

dycosh~2''(?/) . (4.23) 
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In Fig. 6, the amplitude of the bound state is plotted as a function of the brane thickness for each coupling, with 
fi = H . Interestingly, the resulting amplitude is almost independent of the brane thickness a and still finite in the 
thin wall limit. 

Thus, as expected, the divergence of the total amplitude in the thin wall limit arises solely from the KK contribution. 
Regardless, for finite values of cr ~ 0.1 the total amplitude settles down to finite positive values. The result shows 
that the surface divergence for the KK modes can be regularized by introducing a finite brane thickness. This is one 
of the main results of this article. 
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FIG. 5: The total amplitude is shown as a function of 
the brane thickness, a, in the case of d = 2. The thick, 
thick-dashed and dashed curves correspond to the cases 
of ^ = 0, 1/32, 1/16, respectively. 



FIG. 6: The amplitude of the bound state is shown as a 
function of the brane thickness, a, in the case of d = 2, 
with fi = H. The thick, thick-dashed and dashed curves 
correspond to the cases of ^ = 0, 1/32, 1/16, respectively. 
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FIG. 7: The running of the bound state is shown as a func- 
tion of the renormalization scale In fj, in the case of d = 2, 
with o = 0.01. The vertical and horizontal axes show the 
bound state amplitude and log^Q(/i/J/), respectively. The 
thick, thick-dashed and dashed curves correspond to the 
cases of ^ = 0, 1/32, 1/16, respectively. 
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FIG. 8: The relative amplitude of the KK modes to the 
bound state mode is shown as a function of the brane thick- 
ness, CT, for minimal coupling, ^ = 0, for d = 2, with fi — H. 



The bound state amplitude depends on the choice of renormalization scale, ^. In Fig. 7, the running of the scale is 
shown as a function of fi. It is essentially proportional to In /i. The tilt becomes steeper for smaller coupling parameter 
^. There are several possible choices for the renormalization scale, for example, one could choose the expansion rate of 
the brane fi — H or another choice is the brane thickness fi — H/a. We still have no signature from braneworld today 
and therefore no quantity that we can renormalize into. The renormalization scale fi should be determined by future 
observations and/or experiments. In this article we just plot the running of the scale and take the optimal choice 
fi = H iox cases where one has to make a choice. Note that from Eq. (I4.17() the KK amplitude is also proportional to 
ln(/u) with negative tilts. 
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It is also interesting to compare the relative amplitude of the KK modes to the bound state mode. The relative 
amplitude is given by 



(x^(0))kk 
(X^(0)>b. 



KK 



bs 



1, 



(4.24) 



where in the final step we used Eq. (|4.17|) . The result depends on the choice of the renormalization scale /i and brane 
thickness, a. It is meaningful to show the plot for physically reasonable cases. As an example, in Fig. 8, we have 
plotted r as a function of a for the minimally coupled case, ^ = 0, i.e., for tensor perturbations, with fi = H. 

V. KALUZA-KLEIN AMPLITUDE: d = 4 CASE 

In this section, we perform the calculation for the more realistic case of d = 4. The calculation follows in an identical 
manner to the d = 2 case, if only for more tedium. 

A. Amplitude of the KK modes 



In this case the degeneracy factor for the four-sphere {d = 4) is 



dj = 3 (^J' + 2 j + 1) + 2) 
and hence, the zeta function for the total amplitude can be reduced to 

r(n ] ^ - / 0) V U + 3/2) (j + i)(j + 2) 

' ^ 3^ Jcini F{u) ^ [u2 + (j + l/2)2]^772« 



^sinW.-l)|P 



a[/G(e"/2[/^0) - (j + l)(j+2)(j + |) 



(5.1) 



(5.2) 



where we used the properties of bulk propagator Eq. 14.211 . Again, Roman "P" represents taking the Cauchy principal 
value to deal with the pole a,t U — vja. As for the d = 2 case, we divide the total zeta function into a UV piece, 
i.e., for U > 2i//a; and an IR piece, i.e., for Q < U < 2vla. Similarly, the choice of the division is just for later 
convenience. 

To begin with, for the UV piece we shall use the asymptotic formula [23| 



2^ (j + 3/2)(j + l)(j + 2) 
3/2)2]^ 



0- + 3/2)^ 



1 



-T- 



(j + 3/2) 



3 U - + 3/2)2] ^ 6 _ + 3/2)2 

1 



(-ir - 



12(s-l)(s-2)(s-3) 



diB{-U\a,s-2) + 



12(s- 1 



-aa0(-t/2,a,5) 



where 



1 



2(s-l)(,s-2) 
1 



aa0(-[/^a,s- 1) 



a=3/2 



i2r(s) 



r(j + .5 - 1) 



t/-2-'a,Cif(-2j,a) 



a=3/2 



^-2(.-3)g^ r(, + p3) ^_,,^3^^^_^^.^^, 

oo 
J=0 



a=3/2 



a=3/2 



(5.3) 



-"o [U + a)2 + 



(5.4) 
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Then, for the IR piece, we use the following binomial expansions |27l l28|: 



1 g 0- + l)0- + 2)a + 3/2) ^ 1 g rX^^.. ( + 2J - 3, r) - .C.(2. + 2J - 1, 



J!r(s) 



-) 

2' 4 



J = U LW ■ - J ,7=0 

which is valid for the range < U < 3/2; while for the range 3/2 < U < 5/2 we must use 

lf^ (j + l)(j+2)(j + 3/2) _ 1 
[0- + 3/2)2 



(5.5) 



^ ^CHi2s + 2J-3,h-]cH{2s + 2J-l,l[ 



J!r(s) 

j=o ^ ^ 

and finally, for the range 5/2 < U < ^ we have 

1 V (j + l)a + 2)(j+3/2) _ 1 5 

3^0 [(j + 3/2)2 -C/2]^ (9_[;2)« + (^_[^2)^ 



(5.6) 



^ r(s_K/)j^2j 



J=0 



Jiris) 



Ch{2s + 2J - 3, ^) - icH(2s + 2J - 1, 1 



-2s-2,7 



-2s-2J 



,(5.7) 



The total zeta function is obtained from Eq. H4.6|l . 

First, let us focus on the analytic continuation of the UV piece. Some simple manipulations lead to the following 
expression 



Cuv(0,.s) 



1 m2(s-i) 

12 i/2s 



1 



I C/-2«+3 sin[^(s - 2)] [ -daCHiO, a) 

I Is — 1 



a=3/2. 



C/-2«+^sin[7r(s-4)] 



(s-l)(s-2)(s-3) 



d^CH{0,a)\ , + 



U' 



a=3/2 (s-l)(s-2) 



5aCff(-2,a) 



o=3/2 



sin[^(s-3)] 
~' r(j+.s-2)^^_ 



a=3/2 ^ j!r(s) 



a=3/2 



1 [/-2 
(s— l)(s — 2) a=3/2 S— 1 



a=3/2 



J=2 



j!r(s) 



C/"'^9aCH(-2j, a) 1} 

a=3/2J J 



(5.8) 



Like for the d = 2 case, after analytic continuation to s ^ 1, this leads to 



i?^limCuv(0,s) = -ttct 

s— >1 J 
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1v I a 



i^.-^'^(O) 
P.-^-(O) 

p-t^-(o) 



2€-3 



2vla Vp-^-'(O) 



(5.9) 



where u'i(^, a) are the coefficients of the asymptotic expansion in Eq. 14.9|l given by Eq. H4.10|) . for d = 4. 
The IR piece is already finite for s — > 1, and some calculation shows that 



Cir(0,s) 



2^2(^-1) 

37riJ2s(_i). 



■ sin[7r(s - 1)] P 



j7G(e"/2[/^0) 
F(e-/2C/) 
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E 

J=0 

2^1 



j\ns) 



Cni^s + 2J - 3, ^) - ^Ch{2s + 2 J - 1, ^) 



2(s-l) 



7riJ2-(-l) 

CXD 

+ E 

J=0 

OO 

+ E 



— sin[^(s-l)]P / 



^/2^^(7[/G(e'^''/2[/,0) 



F(eW2[/) 



J!r(s) 



i UHi2s + 2J - 3, ^) - JCff + 2J - 1, 1 



sin[7r(s - 1)] P 



5/2 



j=o 



r(£+J) 2,7 
J!r(s) 



F(eW2t;) \{l-mY {f-U^Y 

-2s-2J 



1 



- Ch{2s + 2J - 3, -) - -Ch(25 + 2J - 1, - 



(5/2) 



2s+2J 



(5.10) 



Note that the number of terms depends on the range of U. For 3/2 < 21^ /a < 5/2 the third term should be dropped; 
while both the second and third terms should be dropped if 2z^/(t < 3/2. 
In the s — > 1 limit, as before, just terms with leading order 



CHi2s - 1, ^) - 2(1^ - ^(3/2) + 0{s - 1) , 
contribute to the resulting IR amplitude. Thus, taking the limit s — > 1, we obtain the IR amplitude as 



(5.11) 



H2liniCiR(0,s) = -iaP 



2''/"^^C/3G(e"/2[/,0) 



2l//(T 



dU 



1 



2l//(T 



12 / 



1 



i^(e"/2/7) 

I 2^2 

I r(- -L J 



1 

— crP 
12 



'"/"^^C/G(e"*/2C/,0) 



i)r(| + ^ 



F(eW2[/) 

1 /^/^^3 r(ty+ i)' 



U 



2 , 2 , i)r(-| + ^ + i) \<jI n^ + i) 
r(-f + ^^ + i)r(| + ^^ + i) 1 fv\T{u + \) 



r(| + ^ + i)r(-f + ^^ + i) ^\a)T{v + \ 



(5.12) 



where in the final step Eq. H4.15(l was used. 

Finally, we obtain the total regularized amplitude from Eq. H4.16|l . Furthermore, the KK amplitude is obtained 
by subtracting the bound state amplitude (evaluated in Appendix A) obtained from Eq. (|4.17|) . Note that the KK 
amplitude again has a dependence on the renormalization scale /7. 



B. Results of numerical calculations 



In Fig. 9, a numerical plot of (x^(0))tot is shown. Again, the divergence for the thin wall limit can be seen. The 
power of the divergence can be estimated as follows: the dominant contribution in the thin wall limit comes from the 
first term on the right hand side of Eq. (|5.9() . By changing variables to a; = Ua and following the same steps as for 
d = 2, we obtain 




-1f;Jl-i:«*e.--')^i:q^(^)'-") 

k (9 (^(f ) - + + ^0 - ^^^) + 0(.-) . (5.13) 

In this case, the contribution from the lower bound of the integration does not contribute to any power of a. Thus, 
in the thin wall limit the regularized amplitude diverges as tr^^. This is more divergent than the case of d = 2 and is 
related to the fact that in higher dimensions we need higher powers of UV subtraction. 
The amplitude of the bound state is calculated in Appendix A and is 



^'(x'(o)>bs ^ '^ycosh- 



2u 
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FIG. 9: The total amplitude is shown as a function of the 
brane thickness, a, in the case of d = 4. The vertical axis 
shows logj^Q(— _ff^(x^(0))KK). The thick, thick-dashed and 
dashed curves correspond to the cases of ^ = 0, 3/32, 3/20, 
respectively. 
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FIG. 10: The amplitude of the bound state mode is shown 
as a function of the brane thickness, a, in the case of d = 4, 
with n = H. The thick, thick-dashed and dashed curves 
correspond to the cases of ^ = 0, 3/64, 3/32, respectively. 



,7=2 ^ ^ J=0 ^ ^ 

This is plotted in Fig. 10 and we see that the bound state is almost independent of the brane thickness and still finite 
in the thin wall limit. Thus, like for d = 2, the divergence in the total amplitude arises solely from the KK modes. 

Again, the amplitude depends on the choice of the renormalization scale fi. As we stated in the previous section, 
we have no signature from braneworld as yet and no way to determine the renormalization scale. In certain cases the 
amplitude of the bound state can become negative. In Fig. 11, the running of the bound state amplitude is shown 
as a function of fi. It is basically the same as the case of d = 2; however, a new feature is that negative tilts of the 
running are realized for larger values of coupling ^ which satisfy 

as can be seen from Eq. H5.14I) . The critical coupling parameter in Eq. (|5.15|) is smaller than conformal coupling, 
= 3/16, for any choice of brane thickness, a. This fact means that there always exist coupling parameters which 
realize negative tilts of the running. 

The relative amplitude of the KK to bound state ratio, defined by Eq. (|4.24|) . depends on the choice of renormal- 
ization scale /i. Again as one of the possible physical choices, in Fig. 12, we plot the relative amplitude in the case of 
fi ^ H for minimal coupling, ^ = 0. 



VI. SUMMARY AND DISCUSSION 



We discussed the quantum fluctuations in a thick brane model in order to examine whether or not a finite brane 
thickness can act as a natural cut-off of for the Kaluza-Klein (KK) mode spectrum. The thick brane model we 
examined was supported by a scalar field with an axion-type potential. The thin brane limit of this model is smoothly 
matched to the system of a de Sitter (dS) brane in a Minkowski bulk. As we showed for general d + 1 dimensions, 
this model is classically stable both against classical tensor and scalar metric perturbations (see Appendix B). 

Next, we introduced a test quantized scalar field, into this model and calculated its amplitude. This scalar field 
is assumed to have a zero bulk mass and non-zero coupling to the background domain wall geometry. For simplicity 
we ignore any explicit coupling of x to the geometry; namely, to the geometry or the supporting scalar field (f). A 
particularly interesting case is that for minimal coupling, ^ = 0, which (for d = 4) is equivalent to that of the tensor 
perturbations for a possible low energy realization of general relativity. 
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FIG. 11: The running of the amphtude of the bound state 
is shown as a function of the brane thickness, a, in the 
case of d = 4, with a = 0.01. The vertical and horizon- 
tal axes show the bound state amplitude and \og^Q{^/H), 
respectively. The thick, thick-dashed and dashed curves cor- 
respond to the cases of ^ = 0, 3/32, 3/20, respectively. 
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FIG. 12: The relative amplitude of the KK modes to 
the bound state mode is shown as a function of the brane 
thickness, a, for the minimal coupling, ^ = 0, for d = 4, 
with fj. = H. The vertical axis shows log^Q |r(/i)| given by 
Eq. JTM . 



The squared amplitude of the KK modes was evaluated by using the dimensional reduction approach developed 
in Ref. In order to obtain a well-posed quantum field theory, we introduced a regulator boundary into the set 
up implying that the KK modes become discrete. We worked in Euclidean space rather than the original Lorentzian 
space. For the purpose of the calculation of the amplitude, we used the local zeta function method, where by ^^locaF 
we mean that the quantity is integrated out over the volume of the d-sphere, S'^ (which is d-dimensional dS spacetime 
in Euclidean space), but not over the extra-dimension, z. This quantity is described by a summation over all the KK 
modes and internal modes associated with dS space. Then, given the residue theorem, we can convert the summation 
of all the KK modes into a contour integral representation. The KK modes correspond to poles on the real axis and 
the contour C is taken to enclose them, which is depicted in Fig. 3. Furthermore, as in Fig. 3, it can be deformed to 
the contour C . Then, by keeping the index of the zeta function s large, we can obtain an integral along the imaginary 
axis with two clockwise semi-circles to avoid the bound state poles. Such a contour, C", corresponds to the subtraction 
of the bound state contribution from the total amplitude, i.e., the KK contribution. However, for technical reasons 
it was convenient to integrate along the contour C, depicted in Fig. 4, by adding the bound state contour onto the 
contour C . The quantity obtained from the integration along C is the total amplitude. Thus, after subtracting the 
bound state part, evaluated in Appendix A, we were able to get the desired KK amplitude. 

The bulk propagator can then be decomposed into two parts, where one is regulator-brane independent and the 
other is regulator dependent. Then, by sending the regulator-brane away from the domain wall to infinity, the 
regulator dependent part vanishes and we can take a well-defined single brane limit. 

We decomposed the total zeta function into a high frequency (UV) piece and a low frequency (IR) piece to accomplish 
a successful regularization. Then, we summed these pieces up to get the total amplitude. 

As an exercise, we first calculated the quantum fluctuations for the d = 2 case. For extremely small thicknesses 
the UV piece dominates and exhibits a logarithmic divergence in the thin wall limit. For larger thicknesses the 
total amplitude settles down to finite positive values. Furthermore, we also calculated the amplitude of the bound 
state, which depends on the renormalization scale. However, for a fixed scale we find that its amplitude is almost 
independent of the brane thickness and is finite in the thin wall limit. Thus, only the KK modes lead to surface 
divergences in this limit. In other words, the KK amplitude is regularized by the presence of a finite brane thickness. 
This is one of the main results in this article. We also discussed the running of the bound state amplitude and (of 
particular physical significance) the relative amplitude of the KK modes to the bound state contribution. 

Then, we calculated the same quantities in the more realistic case of d = 4. The main results are very similar to 
the d = 2 case and need not be recapitulated. The main difference is the divergent behavior in the thin wall limit. 
We showed that the KK amplitude diverges quadratically, as opposed to logarithmically for d = 2. We also showed 
that the KK amplitude has an overall negative magnitude.^ Regardless, for d = 4 the amplitude of the KK modes 



This was for our particular choice of renormalization scale, fi = H; however, in general it has an overall negative magnitude unless one 
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is also free of surface divergences for a finite brane thickness. We have no signature from braneworld today and the 
renormalization scale should be determined by future observations and/or experiments. 

In this article we investigated the quantum fluctuations in a particular model of thick braneworld. However, the 
qualitative behaviour of the quantum fluctuations should be independent of the choice of the model. This can be 
understood as follows. For a thick branc model, for which the spacetime is smooth everywhere, there will be no 
divergence. Now if we look at the behaviour of the background solution, Eq. (|2.6(l . when a is sufficiently small, we 
find ^ (j)f){Hz/a) for Hz <^ cr. This is a very general behaviour that one finds at the center of any domain wall 
solution, independent of the global features of the bulk potential. Thus the divergence in the thin wall limit is due to 
the spacetime singularity caused by the divergence of d(j)/dz = at z = 0, which is common to any thick brane model 
supported by a bulk scalar field. 

We are also interested in the amplitude of the bulk tensor metric perturbations, i.e., gravitons. As we mentioned in 
Sec. Ill A the wave equation for the massless, minimally coupled, test scalar field discussed in this article, is equivalent 
to that of the tensor perturbations and we therefore expect a similar result; though an explicit demonstration is left 
for future work. 

In summary, in this article we have shown that a finite brane thickness acts as a natural cut-off for the KK spectrum. 
This fact implies that brane models which have a finite thickness are more plausible than infinitesimally thin ones. 

Some issues remain. In this article, we focused on the fluctuation amplitude just at the center of the wall z — 0, 
i.e., the brane position in the thin wall limit, but the configuration of the fluctuations in the bulk, i.e., the z dependence 
of the amplitude, especially close to the brane, is also important. Moreover, we should also evaluate the Hamiltonian 
density, which is an important quantity in its own right, particularly concerning the back-reaction of the KK modes 
|es| )ecially for small brane thicknesses). Employing the method developed here, it is possible to evaluate such quantities 

The more realistic case of a thick brane model which is embedded in an asymptotically AdS spacetime^ or into a 
bulk with higher dimensions,^ is also of interest. We hope to report on such topics in the near future. 
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APPENDIX A: QUANTUM FLUCTUATIONS OF THE BOUND STATE MODE 

In the following two sub-sections we evaluate the amplitude for the bound state zero mode. The integration here is 
done along the closed contour with the dotted fine as depicted in Fig. 3. 

1. The two-sphere 

First, we note that the bound state for the minimally coupled case is given by 

90 = - , (Al) 
a 

where the bound state zeta function is defined as 

CM = 4/o2(0)L_ + 1/2) (j + 1/2)2 

= 4/o'(0)^^^Cb.(s) . (A2) 



chooses fj, to be very large. 

® In the high energy Umit, Hi 2> 1 is the AdS curvature radius), the bulk effectively becomes Minkowski and hence should reduce to 

the model presented here. 

^ For example, the quantum effects of thin branes with higher spatial dimensions and alike were recently discussed in I3lll32ll33l . 
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Here, fo{z) is the normalized mode function of the bound state. Quite clearly we have a zero mode (by zero mode we 
mean that the lowest eigenvalue Aq is a null eigenvalue, i.e., Aq — 0) and in such a case we have to project out this 
mode to evaluate the bound state contribution. However, in general the bound state varies from the top of the mass 
gap at [v/cr)'^ = 1/4 down io v/a — 0, which is for the massless conformally coupled case. In the following we shall 
focus on a general bound state mass u /a taking care when dealing with the bound state zero mode. 

It is straightforward to evaluate the above ^-function by employing the binomial expansion method which follows 
identically to that of Allen 'i^, see for the case when null eigenvalues are present. Thus, subtracting out the null 
eigenvalue we obtain the following: 



Cbs(s) 



with Cc(s) for S*^ defined by 



1 ^ r(s + J) 

2 ^ J! r(s) 

J=0 ^ ' 



Us) 



©"Cc(. + J)-^,,o(i 



2CH(2s-l,i 



-2s . 



(A3) 



(A4) 



which is the zeta function for the bound state mode in the conformally coupled case, evaluated explicitly in 8] in the 
case of d + 1 = 3. Essentially the minimally coupled case requires summing from J = 1 instead of J = in l)A3p . i.e. 
we have to subtract out the null ei gen value. 

Similar to the case discussed in [23 (section 11.3, Eq. (11.73), pp. 80) there is a pole in the above Hurwitz zeta 
function at s = 1, which can be simply inferred from the relation Eq. (|4.13() . As discussed in ral a suitable way to 
deal with the the pole at s = 1 is to apply the improved zeta function method, described in [35ll36j|. which leads to a 
expression for the amplitude 



(X'(a:))bs = lim-[(s-l)Cbs(s)(x)] . 
s^i as 

Note, the above expression agrees with the usual definition when there is no pole at s = 1. 
Applying the above equation to our case we obtain 



(A5) 



H^{x\Q)),s = 2/2(0)(21n(|)-2V^(l/2)-<5^,oQ) % ^ [2 (^) C/f (2 J + 1, i) - -5^ ,0 Q 



J=i 



Next, we determine /o(0). The normalized bound state solution is 

/2(z) = -Lcosh-2''(x) (^^"dycosh-^'^(y) 



Thus, 



/o(0) = ^ (^"rfycosh- 



'{y) 



(A6) 



(A7) 



(A8) 



Note that for the conformally coupled case — S^ct v ^ Q and therefore, the amplitude of the bound state vanishes. 
This agrees with the result found in 8], for the thin brane case. In fact numerical plots of the amplitude for the 
bound state mode versus the brane thickness show that the bound state mode is independent of the brane thickness. 
Finally, for the bound state mode, we obtain the normalized amplitude as 



i?'(r(o))bs = 



2 In 



-2V;(l/2)-5^^o 



,7=1 



-6, 



(A9) 



This can now be compared with the result for that of the KK modes. 
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2. The four-sphere 



The zeta function for the bound state can be written as 



1 / N 2(s-l) 



where 



. f . _lf; (j+3/2)(j + l)(j+2) 
3^4^J(j+3/2)2-(^/a)2]« 



(AlO) 



(All) 



is the zeta function for a massive scalar field on S''*. For the 5^ geometry, the zeta function for a massless, conformally 
coupled scalar field is given by 



CcW-^[C//(25-3,^)-iCH(25-l,^) 



(A12) 



Thus, the dS zeta function for a general mass can be written as a summation over the massless conformal zeta 
functions (by employing the binomial expansion) 



- , , ^ r(s + J) r /i/\2J ,^ . /3 



J! r(s) 

r(s + J) 
^! r(.) 



J=0 



E 



J=0 



-2s . 



1 /i^\2./ 

.3 

3, 



{Cff (2s - 3 + 2J, ^) - iCff (2s - 1 + 2J, ^)} - ,5^,0 



i[C.(2.-3,^) + .(^)^C.(2.-l,^)-lc.(2.^1,^) 4- 

r(s + J) ^i^^2J^ ^ ^^3, 1 ^ r(s + J) ^Z^N^2J 

£^^7!T(^ 



1 /i^\2 



©^■'C.(2.-3.24)4E 



^! r(.) Va 



Ch(2s + 1,-) 



Cif(2s-1 + 2J,-) 



E 

j=o 



r(s + J) 
J!r(s) 



-2s 



(A13) 



Now, we can evaluate the squared amplitude of the bound state from Eq. ljA5p . The normalization of the bulk 
mode is the same as the d = 2 case and at z = we obtain Eq. IjASp with for d — 4, 

(AM) 



^^-^(Vl + (3-16C)(3fT + 2a2)-l 
The resultant bound state amplitude is 

H'(x'(0))bs = ^^°°dycosh-^^(y) 



In(^) +?Ch(-1,|) + 7^(3/2) -(- 



^ E © - J) - ^)) - E G) 1 

,7=2 ^ ^ J=0 ^ ^ 



1 2 , , , 1 ,3, 
3 + 3^(3/2) + -C.(3,- 



(A15) 



APPENDIX B: CLASSICAL STABILITY AGAINST TENSOR AND SCALAR PERTURBATIONS 



Here, we briefly discuss the stability of the thick brane model both against tensor and scalar perturbations for 
general d-dimensions. 
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1. Tensor perturbations 

We first discuss tlie tensor perturbations about the domain wall background. Here we shall assume a Randall- 
Sundrum (RS) type gauge 1] in which the components of the extra-dimension are zero, i.e., 

ds^ = b^iz) [dz^ + (7^^ + hf,^) dx^'dx''] , (Bl) 

where /i^jy satisfies the usual transverse-traceless gauge about the background dS metric; D^^h^i, ~ /i^^ — 0, where 
C is the covariant derivative associated with 7^,^. 

In this case, the perturbation is separable and we obtain the equation of motion in the bulk direction, which can 
be written in the standard quantum mechanical form as 



dz 



^+Vt{z) iAz)^m'^{z), (B2) 



where °c b{z) '^^^'^h^y and 



{d-lf. d-1, d-l + f 



Vt(z)^^ '-H^ 5 2—. (B3) 

The thin wall limit can be obtained from the limit — > 0, which leads to a system composed of a thin dS brane 
embedded in a flat Minkowski bulk (for quantum fluctuations in such a model, see e.g., ,^). The potential for the 
tensor perturbations in the thin wall limit is then 

Vt{z) = _ (d _ i)H5{z) , (B4) 

where we used for a —t 

^j-—^5{x)- (B5) 

2crcosh [x/a) 

In this limit the solution for the tensor perturbations reduces to the standard exponential form. 

The general solution can be decomposed into a zero mode with mass to = 0, which may realize four-dimensional 
gravity on the brane, and a continuous spectrum of Kaluza-Klein (KK) modes with to > 3/2 (in the five-dimcnsional 
case). Thus, the model is classically stable against the tensor perturbations. 

2. Scalar perturbations 

Next, we discuss the stability of the model against scalar perturbations. We consider a scalar metric perturbation 
of the form 

ds^ = b{zf \{l + 2A\dz'^ + 2Df,Bdx''dz + ("/^,Jl + 27^) + 2D^D^E]dx''dx''] (B6) 



and also a perturbation of the field (j){z) — > (/)(z) + S(t){x), which supports the domain wall. 

In the bulk longitudinal gauge, B = E = 0, the perturbed Einstein equations can be written as follows: 

b' / \ 2 3V 

d{d~i)-n' + (d-i)an + d(d-i)H^n~d(d-i)(-] a = cfy'Scf)' ~ Acp'^ -b^—i 

b \ b / ad) 



-{d-l)Djn' --A) ^ (j)'D^54,, 



(d-i)n" + (d~iy^n' + (d-2)an + (d~2)(d~i)H'^n-(d-i)^A' ~ 2{d~xy-A + uA 

b b b 

= A^^ -4^^' -b'^-^l 



D''Dp{{d-2)n + A] = 0. (B7) 
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The perturbed equation of motion of the scalar field is found to be 

(50" + (rf-l)^(5(^' + n5(/)-2^(0" + (d-l)^<^') + i^-A' + dn'^(^' - b^^6(j) ^ . (B8) 
Next, we derive the evolution equation for the curvature perturbation TZ. By defining 

'2 X -1/2 



*=7^(^j (B9) 

the equation for the curvature perturbations can be reduced to the form 

- + Vs"^ = □d^' (BIO) 

with potential 

. i±^i^y _ 4£ . 

For the dS thick brane case, which is considered in this article, we obtain 

Vs = — -|2[2 + (3(i- 7)ct- (4d-4)cr2] + U + 4(d - 3)cr + (rf^ - I0d + 9)a^] sinh^{Pz)] . (B12) 

4 cosh {/3z) I- L J 

Thus, it is simple to see that at least both for the cases of interest, d = 2 and = 4, Vs > and therefore, the model 
is always stable against scalar perturbations. The d ~ 4 case was originally derived in [20|. 
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